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ABSTRACT 


A  detailed  exp  oration  of  free  electron  laser  (FEL)  theoiy  has  been  done  in  two 
areas.  An  exact  solution  to  the  phase-space  trqectories  in  a  linearly-polarized 
undulator  has  been  obtained  using  a  numerical  simulation.  The  complicated 
phase-space  motion  caused  by  transverse  undulator  deflections  makes  a  rigorous 
derivation  for  trsgectories  difficult,  if  not  impossible.  The  numerical  solution 
extends  the  understanding  of  electron  trsyectories  by  quantitatively  describing  the 
fast  and  slow  components  of  motion.  The  Bessel  ffinction  coupling  coefficient, 
c/o(^)-(7i(^),  describing  the  slow  evolution  is  found  to  be  valid  over  a  broad  range  of 
parameters  even  though  its  derivation  is  approximate. 

A  second  program  has  been  developed  that  provides  a  simple,  quick  diagnostic 
for  accelerator  designers  to  evaluate  how  well  a  simulated  beam  design  will 
perform  as  an  FEL.  The  effect  of  beam  quality  conditions  like  energy,  angular,  and 
positional  spread  are  shown  to  depend  only  on  the  initial  conditions  of  the  beam  at 
the  entrance  to  the  undulator.  This  program  takes  the  six  phase-space  coordinates 
of  the  beam  directly  firom  an  accelerator  simulation  code,  like  PARMELA,  and 
predicts  its  performance  in  an  FEL  system.  This  method  substitutes  for  more 
lengthy,  complex  integrated  simxilations,  like  INEX,  that  require  a  CRAY 
computer. 
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I.  INTRODUCTION 


Since  its  conception  in  1971,  the  free  electron  laser  has  been  the  subject  of 
extensive  theoretical  and  experimental  work  [Refs.  1-3].  While  progress  is  being 
made  on  advanced  designs,  the  fundamental  theory  of  FELs  is  constantly  being 
improved.  There  still  remains  many  aspects  of  the  basic  theory  which  have  not 
been  solved  exactly.  The  essential  physics  of  the  PEL  is  determined  by  the  single¬ 
particle  trajectories  in  the  undulator  and  the  evolution  of  the  optical  fields  [Refs. 
4,6].  Although  FELs  using  linearly  or  circularly-polarized  undulators  work  on  the 
same  principle,  their  electron  phase-space  tr^ectories  are  drastically  different. 
The  exact  trsgectories  in  a  linearly-polarized  undulator  are  complicated  by  motion 
coupled  to  the  transverse  imdulator  deflection  which  causes  the  longitudinal  com¬ 
ponent  of  velocity  to  oscillate  twice  each  undulator  period.  This  causes  difficulty  in 
analysis,  because  these  less-interesting,  fast  longitudinal  oscillations  cannot  be 
precisely  separated  fix)m  the  more-interesting,  slow  bunching  motion  responsible 
for  gain.  Historically,  the  fast  motion  has  been  averaged  into  a  coupling  coefficient 
which  modifies  the  equation  of  motion  for  the  slow  evolution  [Refs.  6-8].  A  small 
error  in  the  coupling  coefficient  could  result  in  a  large  difference  in  the  growth  rate 
or  ultimate  power,  especially  in  high-gain  FELs. 

It  is  important  to  consider  the  validity  of  the  coupling  coefficient  used  to 
describe  the  linearly-polarized  undulator.  The  coupUng  coefficient  uses 
A'[Jo(^>-Ji(^)],  instead  of  just  K  for  the  helical  case,  where  Jo  and  Ji  are  Bessel 
functions  of  argument  K'^I2{\+K‘^).  K  =eBXo/27tmc^  is  the  undulator  parame¬ 
ter,  where  e  is  the  magnitude  of  the  electron  charge,  B  is  the  undulator  rms  field 
strength,  Xo  is  the  undtilutor  period,  and  mc^  is  the  rest  energy  of  the  electron.  The 
averaging  method  uses  an  approximation  to  the  fast  motion  before  taking  a  time- 
average  over  each  period  of  the  tindulator.  However,  all  derivations  of  the  coupling 
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coefficient  solve  for  the  approximate  fast  motion  without  optical  fields  [Refs.  7-12]. 
This  method  is  not  completely  rigorous,  since  the  optical  interaction  is  not  included 
self-consistently.  In  order  to  make  a  more  complete  analysis,  we  derive  the  equa¬ 
tion  describing  the  z  motion  without  approximation  or  assumption. 

With  the  improvement  of  computer  processing  capabilities,  the  last  ten  years 
have  seen  the  emergence  of  simulation  as  a  necessary  tool  for  some  experiments. 
Projects  that  could  cost  millions  of  dollars  are  being  evaluated  for  feasibility  by 
simulations.  As  confidence  improves,  these  methods  are  becoming  more  useful  for 
designers.  For  FEL  experiments,  which  are  by  nature,  costly  and  complicated,  the 
savings  realized  by  simulation  could  be  substantial.  Programs  like  PARMELA, 
used  at  Los  Alamos  National  Labs  (LANL),  are  now  used  to  design  electron  beam 
accelerators  before  projects  are  built.  The  ability  to  evaluate  these  designs  as  they 
perform  as  a  system  is  essential.  Integrated  simulations  like  INEX  currently  per¬ 
form  this  role,  but  as  the  accxiracy  improves,  so  does  the  size  of  the  code,  and 
therefore  the  time  and  money  to  run  them.  Another  feature  of  such  programs  is 
the  large  amoimt  of  input  required  and  prodigious  output  produced.  Both  serve  to 
discourage  casual  use.  A  simple,  quick  diagnostic  that  accounts  for  the  essential 
physics  can  fill  the  gap  between  accuracy  and  simplicity,  and  in  turn  serve  as  an 
accessible  tool  for  the  designer. 

The  effects  of  beam  quality  on  the  performance  of  the  FEL  can  be  predicted  by 
the  initial  conditions  of  the  beam  at  the  entrance  to  the  undiilator.  Electron 
accelerator  sim\ilations  must  keep  track  of  the  six  phase-space  coordinates  of  the 
sample  electrons  as  they  travel  through  the  beam-line  components.  The  output  of 
the  accelerator  simulation  can  be  used  as  the  input  for  the  FEL  performance 
evaluator.  The  effects  of  off-axis  undulator  fields,  which  lead  to  betatron  motion 
[Ref.  13],  and  either  single  or  two-plane  focusing  are  included  to  make  the  descrip¬ 
tion  complete  [Ref.  14].  A  simple  equation  that  converts  initial  conditions  into  the 
meaningful  FEL  dimensionless  phase-velocity  is  used  to  produce  distributions 
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which  describe  the  important  parameters  of  PEL  operation.  This  can  also  be  made 
into  a  standard  'orm  for  use  with  other  previously  developed  models  of  beam  qual¬ 
ity  effects. 
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II.  BACKGROUND 


A  FEL  COMPONENTS 

The  complete  FEL  system  consists  of  three  m£uor  sections  including  the 
electron  beam  accelerator  and  transport  section,  the  undulator,  and  the  optical 
section,  as  illustrated  in  Figxire  1.  The  accelerator  section  is  the  most  complex 
section,  and  includes  not  only  beam-line  components,  but  auxiliary  systems  such  as 
RF  power  and  transport,  vacuum  systems,  beam  fo^sing  and  conditioning 
components,  as  well  as  shielding  and  diagnostic  support.  As  an  amplifier,  the  FEL 
magnifies  a  pre-existing  li^t  beam.  In  the  amplifier  configuration,  the  optics  may 
consist  of  a  seed  laser  and  focusing  elements.  As  an  oscillator,  the  FEL  starts  fit)m 
spontaneous  emission,  and  builds  up  optical  power  until  the  losses  equal  the  gain. 
In  the  oscillator  configuration,  the  optics  of  a  resonator  cavity  is  made  up  oi'  two  or 
more  mirrors,  or  may  take  on  a  variety  of  complex  configurations.  This  thesis 
deals  almost  exclusively  with  the  interaction  inside  the  undulator,  and  assumes  the 
electron  beam  is  provided  with  the  necessary  parameters,  and  that  the  optical 
fields  already  exist  in  well-defined  transverse  and  longitudinal  modes. 


Firure  1.  Elements  of  a  typical  FEL. 


The  transverse  deflection  of  the  highly  relativistic  electron  beam  in  the 
alternating,  periodic  fleld  of  the  undulator  leads  to  classical  spontaneous  emission 
over  a  narrow  spectrum  and  into  a  ?.mall  forward  cone  about  the  longitudinal  axis 
[Ref.  16].  Some  of  this  radiation  may  be  trapped  inside  the  optical  resonator,  and 
provide  the  initial  field  for  an  oscillator.  The  interaction  of  the  electrons  with  the 
combined  fields  of  the  light  and  undulator  cause  microscopic  bunching  on  the  scale 
of  an  optical  wavelength.  This  microscopic  bunching  leads  to  coherent  emission 
and  optical  gain  [Ref.  17]. 

Beam  quality  must  be  sufficiently  good  to  preclude  degradation  of  the 
microscopic  hunching  either  through  energy  or  angular  spread.  A  real  energy 
spread  causes  the  beam  to  "stretch  out,"  because  the  individual  electrons  have 
different  transit  times  through  the  undulator.  If  the  energy  spread  is  large  enough 
to  cause  an  axial  spraad  in  adjacent  positions  which  is  greater  than  an  optical 
wavelength,  then  the  gain  mechanism  will  be  adversely  affected.  A  spread  in  the 
injection  angle  has  a  similar  deleterious  effect.  The  angular  spread  translates  into 
a  spread  in  longitudinal  velocity  as  the  individual  electron  velocities  are  projected 
onto  the  vmdulator  axis.  The  effective  mechanism  is  the  same  as  for  real  energy 
spread  and  causes  a  longitudinal  stretching  that  works  against  the  microscopic 
bunching  needed  for  gain. 

The  electron  beam  needed  for  PEL  operation  may  take  on  a  wide  range  of 
parameters  depending  on  the  use  intended.  Electron  soiurces  include  induction 
linacs,  RF  linacs,  electron  storage-rings.  Van  de  Graaffs,  and  microtrons,  to  name  a 
few.  Typically,  the  beam  is  provided  as  a  series  of  RF  pulses  that  vary  in  peak 
current  fix)m  1  A  up  to  several  kA.  A  direct  current  beam  is  less  common  and 
must  operate  at  lower  current  because  of  power  hmitations,  but  may  still  be  as 
high  as  10  A.  The  beam  radius  may  be  from  less  than  1  mm  up  to  1  cm.  The 
beam  energy  ranges  fix)m  1  MeV  up  to  several  GeV.  Beam  quality  requirements 
depend  on  the  specific  FEL  parameters  for  reasons  discussed  in  the  previous 
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paragraph.  A  more  detailed  discussion  of  the  effect  of  beam  quality  to  FEL 
performance  is  the  subject  of  the  last  chapter  of  this  thesis.  Beam  parameters  can 
be  tailored  to  suit  the  desired  laser  output  quantities  such  as  ultimate  wavelength, 
and  optical  power,  but  because  of  the  cost  and  complexity,  beam  parameters  are 
essentially  fixed  after  the  initial  design.  The  beam  parameters  must  be  combined 
with  the  specific  undxilator  and  optics  parameters  to  complete  the  design  of  the 
laser. 

The  optical  section  for  an  oscillator  usually  consists  of  a  resonator  similar  to 
standard  laser  designs.  This  may  be  a  pair  of  spherical  mirrors  aligned  axially, 
which  sustain  a  transverse  mode  structure  such  as  me  fundamental  Gaussian 
mode  [Ref.  18].  A  typical  resonator  length  is  about  5-10  m  for  an  undulator  that 
may  be  1  to  3  m  long.  For  optical  wavelengths  this  represents  a  mode  separation 
of  X^/2S  =  0.001  A  for  representative  FEL  parameters,  where  X  is  the  optical 
wavelength,  and  S  is  the  resonator  length.  The  gain  bandwidth  is  substantially 
greater  than  this  so  it  is  not  necessary  to  consider  longitudinal  boundary 
conditions.  Resonator  lengths  must  be  chosen,  for  a  given  wavelength,  in  order  to 
ensure  that  the  fundamental  transverse  mode  doesn’t  diverge  too  rapidly  inside  the 
undulator  due  to  diffraction.  This  preserves  a  significant  overlap  between  the 
transverse  dimensions  of  the  electron  beam  and  the  optical  mode.  Beam  quality 
must  also  be  constrained  to  prevent  significant  divergence  for  the  same  reason. 
Angular  spread  will  cause  the  electron  beam  to  diverge  rapidly  or  limit  the 
TnininuiTTi  spot  size  obtainable  for  focusing.  Both  effects  limit  the  amount  of 
transverse  overlap  between  the  electron  beam  and  optical  mode. 

The  undulator  section  is  usually  constructed  of  alternating  magnetic  pole 
pieces  designed  to  present  an  intense,  transverse,  periodic  field  to  the  incoming 
highly  relativistic  electron  beam.  A  typical  peak  field  strength  may  be  as  high  as  a 
few  Tesla,  and  have  a  period  of  a  few  cm.  The  relativistic  contraction  of  the 
undidator  period,  as  seen  by  the  near  light-speed  electron  beam,  causes  a  Doppler 
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shift  in  the  radiation,  and  is  the  mechanism  that  leads  to  short  optical  wavelength 
operation.  Transverse  fields  may  be  constructed  in  a  linearly-polarized  structure, 
which  leads  to  linearly-polarized  light  emission.  As  an  alternative,  the  magnetic 
field  may  be  ciroilarly-polarized  and  can  be  made  firom  helical  electromagnetic 
windings  [Ref.  19].  The  helical  windings  have  a  higher  peak  field  strength  at 
shorter  wavelength  laser  operation.  Undulator  construction  is  a  field  in  itself,  and 
has  been  the  subject  of  much  theoretical  as  well  as  experimental  research.  The 
theory  for  the  helical  undulator  has  the  additional  advantage  that  the 
mathematical  analysis  is  simpler,  and  therefore  is  used  as  an  example  in  many 
theoretical  models.  The  important  differences  between  helical  and  linear 
undulators  is  the  msgor  subject  of  Chapter  III. 

B.  ELECTRON  TRAJECTORIES  INSIDE  A  HEUCAL  UNDULATOR 

In  order  to  discuss  the  principle  of  operation  of  the  PEL,  it  is  important  to 
consider  the  individual  electron  tr^'ectories  inside  the  undulator  without  light. 
Typical  electron  beam  sizes  and  currents  are  are  in  the  Compton  regime  such  that 
the  electrostatic  repulsion  between  the  negatively  charged  electrons  have  a 
negligible  effect.  Furthermore,  the  beam  is  focused  externally  before  entering  the 
undulator.  Under  the  influence  of  the  undulator  magnetic  field,  the  electrons  will 
deflect  periodically  from  side  to  side,  or  "wiggle."  The  ideal  form  for  the  on-axis 
helical  undulator  field  is 

cosikgZ ),  sinikoZ ),  0]  ,  (1) 

were  B  is  the  peak  imdulator  field  strength,  kf,  =  2nJXo  Is  the  undulator  wave 
number  associated  with  the  period  “Kg ,  and  z  is  the  distance  along  the  longitudinal 
axis  of  the  undulator.  The  Lorentz  force  law  becomes 

P  = = -e^xS„  ,  and  ■^(ymc^)  =  0  ,  (2) 

at  at 

where  y  =  1/Vl-fP  is  the  Lorentz  factor,  m  is  the  rest  mass  of  the  electron,  e  =  le  I 
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is  the  magnitude  of  the  electron  charge,  ^  =  OVc ,  and  c  is  the  speed  of  light.  Let 
]?  =  ^1  +  where  the  parallel  component  is  along  the  axis  of  the  imdulator,  and 
the  perpendicular  component  refers  to  the  transverse  direction.  Substituting  the 
field  in  (1)  into  (2),  and  separating  the  components  gives 


4-(7P)  =  — ft,B[-8in(A„2),co8(Ao2),0]  ,  and  ^  =  0  .  (3) 

<u  me  at 

We  see  that  the  electron  energy,  "fne^,  and  the  parallel  velocity,  ^y,  are  constants  of 
the  motion.  From  conservation  of  energy,  the  magnitude  of  the  transverse  velocity 
must  also  be  constant.  The  2  motion  can  be  written  immediately  as  2  =  fift , 
where  the  entrance  to  the  imdulator  is  chosen  to  be  at  2  =  0.  Then,  this  can  be 
substituted  into  the  transverse  equation,  which  can  be  integrated  directly.  Tbe 
parallel  component  can  be  foimd  through  the  coupling  equation,  y~^  =  -  pf. 

Assuming  perfect  injection,  the  exact  velocity  components  can  be  written  as 


)  =  yt  cos(p,^oC< ).  sinCpyfeoCf ),  py]  ,  where 


Pil  = 


1+K^ 


1/2 


(4a) 

(4b) 


and  K  =  eB\,/2ianc^  is  the  dimensionless  undulator  parameter  for  field  strength 
B. 


The  maximum  transverse  defiection  can  be  found  by  an  additional  integration 
and  is 


nt)  = 


ZJtPiY 


[  sin(P|^oCt),  -cos(p||A„c« ),  p,(Cfl 


(5) 


For  typical  FEL  parameters,  if  =  1,  X,,  =  5  cm,  and  y  =  100,  this  leads  to  an 
average  deflection  of  r  =  0.1  mm,  which  is  smaller  than  the  typical  beam  radius. 
Even  if  the  deflections  are  too  small  to  be  visualized  on  the  scale  of  the  undulator, 
they  still  are  responsible  for  significant  spontaneous  emission,  and  FEL  gain. 
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The  spontaneous  emission  mechanism  is  essentially  classical,  and  can  be 
described  by  the  Larmor  formula  giving  the  total  power  emitted  from  a  single 
electron  as  P*  =  2e  (in  cgs)  [Ref.  16].  The  transverse  motion  is  already 

found,  and  the  substitution  gives  P^  =  (.8Tir^l3)(.2'^BoCl8n),  where  =  e^lmc^  is  the 
cleissical  electron  radius.  This  radiation  is  narrowband  and  is  emitted  into  a  small 
cone  about  the  forward  axis  of  motion  with  approximate  angular  width  =  y  This 
radiation  can  be  trapped  in  the  optical  resonator  and  used  to  supply  the  starting 
photons  for  oscillator  operation.  The  interaction  of  these  photons  with  the  electron 
beam  in  the  \indulator  is  the  essential  concept  of  PEL  physics.  In  the  amplifier 
configuration,  the  photons  are  supplied  by  the  seed  laser. 

C.  ELECTRON  INTERACTION  WITH  UGHT 
1.  Description  of  Fields 

The  electron  beam  in  the  undulator  with  light  present  will  be  subject  to 
the  combined  undulator  and  optical  fields.  Some  time  after  startup,  there  will  be  a 
classical  optical  field  in  the  sense  that  the  wave’s  magnitude  and  phase  can  be 
measured  with  arbitrary  precision.  The  optical  fields  are  taken  to  be 

Sr  =E[cosikz-<sit-^),-sinikz-<iit-^XO\  ,  and  {6a) 

Sr  =E[  sin(Ar -orf cos(^ -dit 4^),  0]  ,  (6b) 

where  E  is  the  wave  amplitude,  k  =2xA.  is  the  optical  wave  number  for 
wavelength  X,  co  =  ^  is  the  carrier  firequency,  and  ^  is  the  optical  phase.  In  Figure 
2,  the  net  direction  of  the  force,  found  by  the  vector  relationship  (2),  on  an  electron 
depends  on  its  location  within  a  wavelength  of  light.  Half  the  positions  have  a  net 
force  that  retards  their  motion,  and  therefore  lead  to  energy  loss  which  is  taken  up 
in  the  optical  field  as  the  creation  of  photons.  The  remaining  half  of  the  positions 
gain  energy  firom  the  optical  field. 
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INCOMING 

ELECTRONS 


OUTGOING 

ELECTRONS 


Figure  2.  Vector  diagram  ahowing  fields  acting  on  the  electron  [Ref.  20]. 


Since  the  electron  travels  at  slightly  less  than  the  speed  of  light,  a 
niunber  of  optical  wavelengths  will  pass  over  each  electron  as  it  travels  the  full 
length  of  the  imdulator.  When  exactly  one  optical  wavelength  passes  over  the 
electron  in  each  undulator  period  resonance  is  established.  Mathematically,  this 
can  be  expressed  by  the  condition 


o  _  k 


(7) 
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Using  the  solution  for  p||  =  l-(l+K'^)/27®,  for  y»  1,  this  can  be  inverted  to  find  an 
expression  for  the  resonant  optical  wavelength  given  a  specific  electron  energy 


(8) 


Unless  specified  otherwise,  this  equation  will  be  used  to  discuss  the  resultant 
wavelength  for  PEL  operation.  In  practice,  given  a  specific  electron  energy,  there 
are  a  number  of  factors  which  determine  the  ultimate  optical  wavelength,  but  it  is 
roughly  constrained  to  be  in  the  vicinity  of  expression  (8). 

This  relation  also  demonstrates  some  of  the  advantages  of  the  PEL. 
Because  factors  like  the  electron  beam  energy,  the  undulator  field  strength,  and 
the  undulator  period  can  be  modified,  the  operating  wavelength  of  the  PEL  can  be 
tuned  over  a  wide  range.  Also,  because  there  is  no  internal  mediiun  other  than 
electrons  and  co-propagating  lig^t,  there  is  no  physical  hmitation  to  the  power 
density  of  the  laser.  Only  external  components,  especially  the  optics  tend  to  limit 
the  power  of  operation.  These  are  some  of  the  most  desirable  features  of  the  PEL. 


2.  The  Wave  Equation 

Electrons  that  are  close  to  one  optical  wavelength  apart  will  experience 
identical  forces.  Therefore,  a  coordinate  transformation  that  follows  an  electron 
with  resonant  velocity  for  the  chosen  wavelength  is  .  This 

represents  the  phase  of  the  electron  within,  roughly,  one  optical  wavelength.  The 
physical  situation  is  identical  over  a  period  of  =  2x.  A  phase  velocity 
v=L[ik  1  describes  how  much  faster  or  slower  the  electron  is  than  its 

resonance  velocity,  where  v  =  0.  In  order  to  follow  the  growth  of  the  optical  field  it 
is  necessary  to  develop  a  wave  equation.  A  classical  wave  equation  may  be  written 
by  summing  the  individual  electron  currents. 

c 
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where 


□2  =  v2-JL_^ 

c2  af2 


is  the  d'Alembertian  operator,  and  A  is  the  vector  field 


potential,  and  is  the  ciurent  density  [Ref.  16].  An  appropriate  vector  potential  is 
A*  =  (£/^)[  sin(Az-cirf+<|)),  cos(A2-<oi +<])),  0].  Assuming,  the  amplitude  and  phase 
vary  slowly,  and  keeping  only  terms  with  single  derivatives,  the  left  hand  side  can 
be  written  as 


dz  C  dt 


The  total  beam  current  can  be  represented  as  the  s\im  of  the  single-particle 
cinrrents  Jf  = -ec0i5^^^(?^-r5),  where  follows  the  trajectory  of  the  particle. 
Using  the  motion  of  the  electron  in  the  undulator  field,  the  form  for  can  be 
inserted.  Since  the  physical  situation  over  one  wavelength  of  light  in  Figure  2  is 
completely  periodic,  it  is  only  necessary  to  sum  over  a  period  of  X.  The  result  is 
[Ref  17] 

O 

a  =-j<e~*h>  , 

where  the  complex  dimensionless  field  envelope  and  current  are 

.  4KNeKLEe^^  .  _  8N{enKLfp 

®  9  9  t  J  * 

Y^nc  ymc 

E  is  the  peak  field  strength,  and  p  is  the  electron  beam  particle  density.  The  open 
circle  denotes  a  time  derivative  with  respect  to  the  dimensionless  time  variable, 
x  =  ct/L,  and  the  brackets  indicate  a  phase  averaging  over  the  sample  electrons 
within  a  wavelength  of  light.  If  many  electrons  are  uniformly  spread  in 
longitudinal  position,  then  the  phase  average  ,  will  be  zero,  and  there  will 

be  no  optical  field  growth.  In  order  to  understand  the  FEL  interaction,  first  it 
must  be  understood  how  electrons  evolve  firom  their  initial  positions  within  a 
wavelength  of  li^t. 


(9) 


(10) 
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3.  The  Pendulum  Equation 

To  solve  for  the  motion  of  electrons  with  lig^t,  again  write  the  full 
equations  of  motion,  the  complete  Lorentz  force  law  is 

f  =  ^(TfmclO  =  -e[£  +  ( px^)]  ,  and  -^(ymc^)  = 

at  at 

Inserting  the  full  combined  fields  of  the  undulator  and  optical  fields,  the  transverse 
motion  with  perfect  injection  can  be  found  again  by  direct  integration. 

^4.  =  coa(koZ ),  sinikoZ ),  0]  -  8in(to-<of  +<>),  cos(^  -int  -h|>),  0]  ,  (11) 

where  A  -eE  }J2ianc^,  analogous  to  the  definition  of  if .  The  energy  equation  is 


The  form  for  can  be  inserted  to  yield  the  energy  change  equation, 

Y=  (13) 

Using  the  relation  y  an  equation  in  the  longitudinal  variable  fig  can  be 

written.  By  rewriting  the  equation  in  terms  of  the  dimensionless  variable  C,  it  can 
be  shown  that  individual  electrons  follow  motion  described  by  a  pendulum  equation 

^  =  V  =  lo  lcos(C4^)  (14) 

As  electrons  evolve  along  phase-space  orbits,  like  a  pendulum,  they  move  in  phase 
in  such  a  way  as  to  either  drive  or  diminish  the  optical  field  through  (9).  Phases 
which  make  the  argument  of  cosine  between  -n/2  and  iU2  cause  a  gain  in  electron 
energy  and  a  corresponding  loss  to  the  optical  field.  The  other  half  of  the  possible 
phases  lead  to  optical  field  growth.  Using  the  pendulum  equation  it  is  possible  to 
see  by  mathematical  means  what  was  said  earlier  for  a  uniform  distribution  of 
initial  phases,  half  the  electrons  gain  energy  and  half  lose.  But,  imder  the 
influence  of  the  combined  fields,  the  electrons  evolve  and  can  "bunch"  at  desirable 
phasm  and  cause  gain.  The  understanding  of  what  leads  to  hunching  of  electrons 
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on  the  scale  of  an  optical  wavelength  is  the  key  to  understanding  the  FEL.  The 
combination  of  the  pendulum  equation  and  the  wave  equation  form  a  complete 
description  of  the  FEL,  and  can  be  simultaneously  solved  to  show  the  effects  of 
various  initial  conditions  on  operation. 

Figure  3  shows  the  phase-space  evolution  of  electrons  in  a  helical 
undvilator  with  h^t  present,  obtained  by  integrating  (9)  and  (14)  [Ref.  17].  The 
initial  optical  field  is  =  3,  and  the  initial  phase- velocity  is  =  0  (at  resonance). 
The  large  picture  shows  the  phase-space  evolution  of  the  sample  electrons.  The 
dots  change  in  color  firom  light  to  dark  as  they  travel  down  the  undulator.  The 
dimensionless  time  variable  runs  firom  t  =  0— »1.  The  sample  electrons  are  initially 
chosen  to  be  uniformly  distributed  in  C.  The  solid  line  represents  the  separatrix 
which  divides  phase-space  into  closed  and  open  orbits.  The  separatrix  can  be 
foiind  by  an  elementary  analysis  of  the  pendvdum  equation. 

Phase-space  orbits  can  be  fo\md  from  the  pendulum  equation,  (14),  by 
integration.  First,  multiply  both  sides  by  2v  to  make  it  a  perfect  differential,  so 
that 

2w  =  2lo  ICcos(0  , 

where  the  optical  phase  is  neglected.  This  can  be  integrated  to  yield  the  equation 
of  the  phase-space  orbits 

v^-Vo^  =  2  lo  l[sin(C)-sin(Co )]  ,  (15) 

where  v„  and  ^  are  the  initial  phase-space  coordinates.  The  separatrix 
corresponds  to  the  phase-space  path  that  passes  through  ^  =  -x/2,  and  v^  =  0. 
The  resulting  equation  for  the  separatrix  is 

Vggp  =  ±V2la  l(sin(  C  HD  ,  (16) 

which  has  a  peak  value  of  Voj,^^  =  ±2Va7-  Points  which  start  outside  the  region 
enclosed  by  the  separatrix  have  "open”  orbits,  and  points  inside  have  "closed" 
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-11/2  ;  3lt/2  0  t  1 


Figure  S.  Phase-space  evolution  for  monoenergetic  electrons  on  reso¬ 


nance. 


orbits.  There  are  fixed  points  at  ^  =  -)c/2,  3x^2,  and  ic/2.  The  point  at  ^  =  x/2  is  a 
stable  fixed  point  which  can  be  demonstrated  by  considering  small  deviations 
C  s  1C/24X ,  where  x«x.  The  substitution  of  this  into  the  pendulum  equation  yields 

OO 

X  =  la  Ioos(x/24x  )  =  la  I[  cos(ic/2)co6(x  )  -  sin(x/2)sin(x )]  =  -la  lsin(x )  . 

OO 

Using  the  small  angle  approximation,  sin(x)  =  x,  the  equation  becomes  x  =  -la  lx 
which  gives  simple  harmonic  motion  which  is  stable.  The  points  at  -n/2  and  37C/2 
can  be  shown  to  be  unstable  [Ref.  21].  In  Figure  3,  all  points  start  inside  the 
separatrix,  and  therefore  have  closed  orbits.  Also  displayed  are  the  gain  and  phase 
evolution  through  the  undulator.  The  gain  and  phase  are  found  from  solving  the 
wave  equation,  (9). 
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Figure  4  shows  the  effect  of  starting  above  resonance.  The  electrons  are 
started  with  =  2.6,  which  corresponds  to  the  maximum  gain  condition,  which 
will  be  demonstrated  below.  Some  electrons  are  trapped  in  closed  orbits,  and 
others  are  in  open  orbits.  The  open  orbits  move  into  the  next  wavelength,  but 
wrap  aroimd  because  periodic  boundary  conditions  are  used.  The  gain  is 
substantially  higher  and  the  optical  phase  change  is  less,  because  the  electrons 
bunch  near  ^  =  ic. 


Figure  4.  Phase-space  evolution  fw  monoenergetic  electrons  above  reso¬ 
nance. 
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conditions,  and  is  shown  in  Figvire  5.  This  shows  the  maximum  gain  occurs  at 
Vo  =  2.6.  The  gain  spectrum  is  anti-symmetric  about  Vq  =  0,  and  is  characteristic 
of  small  signal  (low  laser  ligd^t  intensity)  and  low  gain  FELs. 


****  Gain  apnctrua  **** 

HbSO 

0.1144 

0.0 

-0.1144 
0.066 

0.0 

-0.066 

10  10 


Figure  5.  Small  signal,  low  gain  spectnun. 
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III.  ELECTRON  DYNAMICS  IN  A  LINEARLY-POLARIZED 

UNDULATOR 

A.  MOTION  WITH  NO  OPTICAL  FIELD 
1.  Equations  of  Motion 

The  previous  discussion  is  based  on  the  helical  undulator  design.  As 
stated  earlier,  there  is  an  advantage  in  the  analysis,  because  the  motion  is  simpler. 
In  this  section,  the  exact  motion  in  the  linearly-polarized  undulator  is  found.  To 
construct  a  systematic  approach,  it  is  best  to  start  with  the  motion  of  electrons  in 
the  presence  of  only  the  undulator  field.  The  optical  field  will  be  considered  later, 
after  the  method  of  approach  has  been  established. 

In  the  helical  case  with  perfect  injection,  the  electrons  travel  down  the 
axis  of  the  undulator  with  a  constant  velocity,  P||.  The  transverse  motion  is 
constant  in  magnitude,  but  changes  direction  smoothly  in  a  circular  fashion.  This 
makes  the  analysis  simple,  because  the  acceleration  is  constant  in  the  equations, 
making  it  possible  to  solve  a  set  of  coupled  differential  equations.  In  the  linear 
case,  the  acceleration  is  not  constant,  but  varies  as  the  electrons  deflect  from  side- 
tO'Side  in  one  transverse  plane.  This  leads  to  difiBculty  in  separating  the  equations 
of  motion.  In  order  to  consider  the  motion,  first  consider  an  ideal  linearly-polarized 
magnetic  field  on-axis 

^  B[  0,  sin(koZ),  C:  ,  (17) 

where  B  is  the  peak  magnetic  field.  Since  the  magnetic  field  varies  sinusoidally 
over  each  period,  the  rms  value  is  E  =  B  /V2.  Assign  cartesian  coordinates  to  the 
undulator  with  the  origin  at  one  end,  the  z  axis  along  the  longitudinal  axis  of  the 
imdulator,  and  the  magnetic  field  in  the  y-plane.  The  electron  motion  will  be  in  the 
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x-z  plane.  The  equations  of  motion  are 


4(TP.)»P.-^sm(»„z)  . 

tu  me 

i<A)  =  o  . 

at  me 


and 


(18a) 

(18b) 

(18c) 

(18d) 


The  x-component  is  a  perfect  derivative,  and  can  be  integrated  directly  to  yield 

P*  = -^^^008(^0-2 )  .  (19) 

where  K  =  eBXof^ione^'  lii  this  case,  perfect  injection  at  z  =  0  corresponds  to  the 
electron’s  x  component  of  velocity  at  a  maximum  in  the  negative  direction.  The 
energy  is  a  constant  of  the  motion,  neglecting  losses  to  spontaneous  emission, 
because  the  magnetic  force  acts  at  right  angles  to  the  electron’s  motion.  Perfect 
ixgection  also  assumes  that  the  other  transverse  component  py  =  0.  Therefore,  the 
perfectly  ii\jected  beam  must  be  sent  into  the  undulator  at  an  angle 
0  =  tan“Hp*/P,)  =  P*/P*-  Since  the  longitudinal  component  P^  =  1,  then  0  = 

For  typical  values  of  FEL  parameters,  K  =  1,  y  =  100,  so  that  0  ~  0.01  rad. 

Substitution  of  the  solution  for  p,  into  the  equation  for  p^  results  in 

ojri 

Px  =  — 2-AoCCOS(Ao2)8in(*o2)  •  (20) 

r 


This  can  be  made  into  a  perfect  differential  by  mxiltiplying  both  sides  by  2^g. 
Then, 


2^cosHkoZ) 

.  r 


,  giving  ^^  =  -2^cosHkoZ)  +  C  (21) 


where  C  is  a  constant  of  integration  that  accoimts  for  initial  conditions.  At  z  =  0, 
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Pi,  =  l-<l+2Ji:2)/y*,  so  that  C  =  I-7  After  some  algebra,  this  can  be  separated 
into  a  form  which  can  be  integrated 


_ dz _ 

VPzo  +  {^'^lf)&vci\koZ ) 


=  cdt 


This  can  be  integrated  and  the  solution  found  in  an  inverted  form,  that  is,  with  t 
as  a  ftmction  of  z ,  as 


t{z)  = 


K 

2iic 


F(9JiI)  , 


(22) 


where  0  =  k^z ,  and  F  is  an  incomplete  elliptic  integral  of  the  first 

kind. 


2.  Phase-Space  Description 

While  this  result  is  pxnely  analytic,  it  is  of  restricted  use,  because  FiBM) 
is  a  tabulated  function.  The  actual  trsyectories  will  be  evaluated  later  by 
niunerical  means.  As  discussed  in  the  helical  case,  the  important  characteristic  of 
the  electron  motion  is  bunching  on  an  optical  scale,  emd  so  the  natural  choice  for 
description  are  the  dimensionless  phase-space  variables,  ^  =  (.k^  +k  )z  -oat ,  and 
V  =L[(A<,+^)P*-A].  Although  there  is  no  light  present,  it  is  still  possible  to  make 
the  discussion  in  terms  of  the  wavelength  defined  by  the  resonance  condition  (7). 
In  the  helical  case,  the  electrons  wo\ild  appear  as  steady  points  in  phase-space, 
(^,v)  that  do  not  evolve.  In  the  linear  case,  the  side-to-side  undulator  deflection 
causes  the  electrons  to  orbit  about  these  steady  points  at  twice  the  undulator 
periodicity.  To  illustrate  this,  first  write  the  complete  expression  for  the 
longitudinal  velocity 

P*(z )  =  VpI,  +  {m'^lf)s,m\Kz )  .  (23) 

In  the  regime  where  {Klyf'cl,  this  can  be  approximated  by  a  binomial  expansion. 
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Pz(2 )  =  P«>  +  ^s\rL\koZ ) 

r 


(24) 


Over  the  full  range  of  possible  arguments  for  z ,  varies  by  Ap^  =  ±K^I2'^.  This 
will  caiise  a  variation  in  phase  velocity  v  by  Av  =  ±L  (A +Ao  lAP^  =  2kNK^I{1+K^). 
For  a  typical  value  of  /ir=l,  and  N=5Q,  this  means  a  variation  of  Av  =  ±5071  during 
each  half  period  of  the  imdulator.  This  is  50  times  greater  than  the  natural  gain 
bandwidth  of  2iil 

The  extent  of  motion  in  ^  cannot  be  computed  without  an  approximate 
solution  of  the  elliptical  equation.  A  perturbation  in  powers  of  (K/y)^  to  first  order 
will  siiffice.  The  result  is 


zit)~^gCt  + - ^cos(2P2^oC<)  +  •  •  >  (25) 

SiC[ 

where  average  value  over  a  half-period.  This  can  be  computed 

analytically  by  the  results  above,  and  is 


p  _jeEL_ 

2F(7C/2,M) 


This  results  in  a  variation  of  Az  =  ±X<,iir^/27iy^. 


)  can  be  computed  directly  as 


(26) 


C(i)  =  (A+Ao)z(i)-(«rf  =^(t>f^cos(2a)oO  ,  (27) 

where  l^it)  =  ik-¥k„)^gCt-wt ^  %=K^I2(.1+K\  and  ©o  =^gkoC.  The  variation  in  ^  is 
A^  =  (/s+Ao)A2  =  ±K'*/2(l+ir*).  For  a  typical  design  with  if  =  1,  A^  =  ±1/4.  For  the 
periodicity  of  2x  in  this  variation  has  a  width  of  7t/6.  This  is  still  only  a  fraction 
of  a  wavelength,  so  it  is  not  expected  to  destroy  the  microscopic  bunching 
mechanism  needed  for  gain.  But,  A^  is  responsible  for  harmonics. 

A  numerical  solution  to  the  phase-space  motion  of  electrons  in  a  linearly- 
polarized  \indulator  is  shown  in  Figure  6.  The  FEL  parameters  used  are  N  =  50, 
X  =  5  cm,  y  =  100,  =  1,  y  =  0  (not  needed  unless  fight  present),  and  Oo  very  small 
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(no  Periodic  boundary  conditions  have  been  used  again.  The  motion  is 

symmetric,  and  the  variations  agree  with  the  approximate  values  computed  above. 

Although  the  sampled  electrons  are  equally  spaced  on  the  average,  in 
Figure  6  they  are  not  actually  equally  spaced  longitudinally  as  they  travel  through 
the  undulator.  The  electrons  are  at  difierent  phases  in  their  elliptical  orbits 
compared  to  adjacent  electrons.  These  differences  are  obscured  by  the  periodic 
nature  of  the  orbits,  however,  it  can  be  observed  by  watching  the  simulation  as  it 
runs.  The  transverse  deflection  slows  them  down,  and  they  tend  to  bunch  as  they 
reach  their  maximum  deflection.  They  spread  out  again  as  they  move  back  to  the 
axis.  Ihis  "artificial"  bunching,  in  the  sense  that  it  does  not  come  from  the  optical 
fields,  requires  some  care  in  selecting  the  electron’s  initial  conditions.  It  does  not 
suffice  to  take  a  uniform  distribution  as  a  starting  point.  The  electrons  must  be 
presented  to  the  undxilator  one  at  a  time,  j\ist  as  it  occiurs  in  the  real  beam.  This 
reqiiires  an  extreme  amount  of  precision  until  all  the  sample  electrons  are  in  place, 
then  the  simulation  can  proceed  as  usual.  Failure  to  provide  for  this  eflect  can 
lead  to  a  non-physical  pre-bunching  which  would  create  false  gain.  Also,  it  was 
found  that  great  care  needed  to  be  taken  in  the  calculations,  because  of  the  fine 
scales  involved,  which  were  on  the  order  of  one  part  in  a  million.  Although  the 
electrons  move  above  and  below  resonance  by  large  amounts,  they  must  return  to 
the  exact  same  conditions  after  each  half-period  of  the  undulator.  In  practice,  to 
calculate  this  requires  double  precision  throughout,  and  a  fourth-order  Runge- 
Kutta-Nystrom  integration  technique  [Ref.  23]. 
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Figure  6.  Exact  phase-space  tr^ectories  of  electrons  in  a  linearly- 
polarized  undulator  with  no  lii^t  present.  The  optical  wavelength  is 
determined  by  the  resonance  condition. 


B.  MOTION  IN  THE  PRESENCE  OF  OPTICAL  FIELDS 

1.  Derivation  of  the  Exact  Equations  of  Motion  in  Ideal  Fieldfg 

Assume  now  that  there  is  an  actual  optical  held  present,  due  either  to  a 
seed  laser  (amplifier),  or  by  captured  spontaneous  emission  (oscillator).  In  the  case 
of  the  oscillator,  with  the  assumption  of  plane-waves,  the  fields  must  have  the 
same  polarization  as  the  motion  of  the  electron.  The  far-field  radiation  will  be  due 
solely  to  the  acceleration  fields,  in  which  case,  the  electric  field  vector  will  be  in  the 
same  plane  as  the  electron’s  motion,  and  the  magnetic  field  at  right  angles  to  it. 
Fields  that  meet  all  of  these  requirements  are 
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~E[  coeOa-<sit-^),  0,  0]  and  Sji  =£[  0,  cos(Acz -<i)f -Ht>),  0]  ,  (28) 

where  E  is  the  electric  field  magnitude,  k  =  2tU\  is  the  optical  waveniunber  for 
wavelength  X,  (D  is  the  carrier  fi:equency,  and  ^  is  the  optical  phase.  The  equations 
of  motion  are 

^(■ymc  Pi,  )  =  -€£cos(fa-0itf  4^) +  cP^[Bsin(feo)  +  ^<»s(fe2-<«)^+^)]  ,  {29a) 

at 

^(•ymcp,)  =  0  ,  (29b) 

-^(‘ymcP*)s=-eP,[Bsin(Ao2)  +  ^cos(fe-<iirf+4>)]  ,  and  (29c) 

at 

■^(ymc )  =  -c  P,B  cos(^  -<at  -hJ))  (29d) 

at 

The  transverse  x-oomponent  equation  can  be  written  as 

4-(‘yP* )  =  ^  cos(v)v  +  >/2K’sin(e)e  ,  (30) 

at 

where  K  ^eBX^I27anc^,  A=eEXl2iimc\  B-koZ,  B  =  koZ,  v*  and 

^zzki-<o.  Note  that  the  time  derivative  of  the  optical  phase  is  ignored,  because 
low  gain  is  assiuned.  (30)  is  a  perfect  derivative,  and  can  be  integrated  directly  to 
give 

P,  =  -V2— cos(0)  +  V2— sin(y)  ,  (31) 

Y  Y 

where  perfect  ipjection  has  been  assiimed.  This  allows  an  exact  expression  for  the 
energy  change  to  be  foxmd  as 

Y  =  ~  cos(0)co8(\(f)  -  A  ^sin(y)cos(  v)j  (32) 

These  two  equations,  (31)  and  (32),  along  with  Y~^  =  1-Pz^-Px  >  are  sufficient  to  find 
a  dififeiential  equation  for  the  z-component.  After  some  rearrangement,  the 
equation  reads 
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P,  (QA^P^-l)sin(2y>fii:A/.(P,  )co8(frfv)-W^(P»  )cos(9-v>HDoiir  W29) 

1-P*  l+2i:2co8^0MiiCAcos(0)sin(v>f2A2sin2(\|/) 

where  <0o  =koC  and  =  (Qo±a)(px-l).  T^thout  A=0,  this  reduces  to  the 
equation  of  motion  vdth  the  undulator  fields  alone,  (20),  and  leads  to  the  elliptic 
functions  previously  derived.  The  complete  equation  as  it  appears  cannot  be  solved 
analytically.  As  an  approximation,  (33)  can  be  shown  to  reduce  to  the  penduliun 
equation  in  the  regime  where  K/y  is  small. 

2.  Pendulum  Equation  for  the  Linearly-Polarized  PEL 

In  the  equation  for  the  energy  change  y,  consider  the  term  compared 
to  KA .  The  value  of  A  depends  on  the  optical  field,  a ,  which  is  limited  by 
saturation.  The  strongest  optical  field  forces  the  electrons  into  about  one 
synchrotron  oscillation  during  their  pass  through  the  undxilator.  Synchrotron 
motion  is  the  periodic  oscillation  in  the  longitudinal  direction  as  compared  to 
betatron  motion  which  is  in  the  transverse  direction. 

The  analysis  of  the  phase-space  motion  shows  that  there  is  a  stable  fixed 
point  at  ^  =  x/2.  The  electrons  in  closed  orbits  follow  nearly  circular  paths  about 
the  stable  fixed  point.  Examine  small  deviations  ^  =  7c/2-4<x ,  where  x  cn: 

OO 

Substitution  of  this  form  into  the  pendulum  equation  yields  x  =  la  lcos(7C/2+x ). 
Use  the  trigonometric  identity  cos(x/24x )  =  cos(x/2)cos(x  )-sin(x/2)sin(x )  =  -x .  This 
gives  a  harmonic  equation  x  =  -la  lx  =  -v^x,  where  v,  =  VloT  is  the  dimensionless 
synchrotron  firequency.  The  assumption  of  approximately  one  S3mchrotron 
oscillation  during  a  pass  through  the  undulator  implies  that  v,  =  2x  =»  la  I  =  47c^. 

If  this  estimate  for  strong  fields  is  substituted  into  the  definition  of  la  I, 

then 

Iq  ^  4v^  .  (34) 

Using  the  resonance  equation  X  =  \oil+K^)/2^,  we  find  that 


A  _  1  a+K^) 
K~  4N^  K'2 


(35) 


For  K=l,  AlK  ~  1I2N^,  so  that  even  for  a  small  number  of  periods,  say,  N:t>10,  the 
ratio  of  terms  is  limited  to  A /K  <  1/200.  For  more  conventional  designs,  the  ratio 
is  even  smaller  yet,  and  therefore,  the  effect  of  neglecting  the  A^  term  is  minimal. 

It  is  now  possible  to  write  the  approximate  energy  equation  as 

7  =  -^Aiircos(e)cos(v)  .  (36) 

By  use  of  a  trigonometric  identity  [Ref.  23],  it  can  be  rewritten  in  a  way  which 
separates  the  fast  and  slowly  evolving  terms 

Y  =  -^AiTfcosCv+e)  +  cos(v-0)]  .  (37) 


The  arguments  can  be  rewritten  as 


V+6  =  (A+A!o)z(i)-<iirf  =  and  ^-0  =  (k+ko)z(t)-iiit-2kot  =^(t)~20  . 

From  (27),  yfO  =  ^(i  )44cos(26),  and  \|^-6  =  X,(t  )+^cos(20)-29.  The  cosine  terms  can 
be  expanded  into  products  of  slowly-varying  and  fast-varying  terms  by  applying 
another  trigonometric  identity  so  that 

cos(v46)  =  cos(^)cos(4cos(20))  -  sin(^)sin(^cos(20)  ,  and 

co8(y-6)  =  cos(^)cos(20-^cos(20))  +  sin(^)sin(20-^cos(20)) 

These  terms  can  be  averaged  over  a  period  of  the  undulator,  assuming  that  Z, 
varies  slowly  and  its  value  remains  essentially  constant  over  one  period.  Over  one 
period,  the  argument  0  varies  from  0  to  2k.  For  an  arbitrary  function  of  0,  the 
average  value  over  one  undulator  period  is 

<»«)>  = . 


The  average  energy  change  rate  then  can  be  evaluated  directly  from  the  following 
integrals  [Ref.  23]: 
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(38a) 


1 

—j^  cos(^cos(20))de  =  Jo(4)  . 

cos(2e-4cos(20))de  =  -T/i(^)  ,  (38b) 

2k 

—  f  sin(^cos(20))d  0  =  0  ,  and  (38c) 

2it 

8in(20-^co8(20))d0  =  O  ,  (38d) 

where  Jf,  and  t/j  are  zeroth  and  first-order  ordinary  Bessel  functions.  Substitution 
of  these  integrals  into  the  averaged  energy  rate  of  change  equation  gives  the 
following  result  [Ref  6] 

<  V  >  =  -AK[Jo  (4)-Ji(^)]cos(C-h1>)  .  (39) 

Y 

This  differs  from  the  helical  case  only  in  the  replacement  of  K  by  iir[Jo(^)-«f 
The  argument  ^  varies  firom  zero  up  to  maximum  value  of  1/2.  The  Bessel  function 
coefficient  therefore  varies  from  a  value  of  1  at  lir=0,  down  to  0.72  as  K-*oo,  This 
can  be  carried  through  everywhere  K  appears,  and  therefore,  it  is  the  only 
modification  necessary  to  derive  the  pendulum  equation  for  the  linearly-polarized 
undulator.  The  slowly-evolving  phase-space  motion  in  the  linearly-polarized 
undiilator  is  then 

OO 

C  =  lo  lcos(C-H|>)  ,  (40) 

where  \a\  =  AnNeK[J„i%)--Ji{i)]LEhfhnc^.  All  of  the  previous  results  from  the 
phase-space  evolution  in  the  helical  undulator  can  be  used  to  describe  the  linearly- 
polarized  imdulator. 

While,  this  result  is  essentially  correct  in  the  regime  where  K/y  is  small, 
this  "standard”  derivation  has  a  number  of  mathematical  problems.  Namely,  a 
number  of  factors  were  assumed  to  be  essentially  constant  over  one  period  of  the 
undulator.  Secondly,  approximate  solutions  to  the  z  motion  are  used  from  the 
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"no-light"  solution  obtained  in  the  earlier  section.  This  cannot  be  complete  because 
the  effect  of  he  light  on  the  motion  is  neglected.  This  is  not  a  rigorous 
perturbation  scneme,  becaxise  the  end  resiilt  contains  all  powers  of  the  optical  field 
in  the  trigonometric  argiunent  in  Other  derivations  proceed  along  the  same 
lines,  and  while  more  sophisticated,  suffer  firom  the  same  inconsistency.  For  this 
reason,  it  is  important  to  return  to  the  exact  differential  equation  (33),  and  seek  an 
alternative  means  of  solution.  A  way  of  accomplishing  this  with  a  minimum  of 
difficulty  is  numerical  simulation.  Since  (33)  cannot  be  solved  rigorously,  not  only 
is  niunerical  simulation  easier,  it  is  more  appropriate  and  can  extend  the  solution 
into  regimes  where  the  JT/y  1  condition  is  not  met. 

C.  NUMERICAL  SOLUTION  TO  THE  EXACT  EQUATIONS  OF  MOTION 
1.  Time-Averaged  Phase-Space  Triuectories 

Equation  (33)  can  be  solved  exactly  by  numerical  means.  This,  however, 
is  not  a  complete  description  of  an  FEL  because  the  field  evolution  must  also  be 
included.  Because  of  this  necessity,  the  solution  shall  be  restricted  to  the  low  gain 
case  such  that  the  optical  field  evolution  can  be  described  by  the  slowly-varying 
wave  equation  (9),  and  where  the  time  derivative  of  the  optical  phase  ({>  can  be 
ignored  in  the  Lorentz  force  equation.  There  are,  however,  no  restrictions  on  the 
values  of  K  and  y.  The  exact  motion  can  be  found  by  integrating  (33),  and 
displaying  the  results  in  the  traditional  phase-space  description  (C,v),  which  was 
developed  earlier. 

Figure  7  illustrates  the  phase-space  motion  of  a  single  electron  in  the 
presence  of  an  optical  field.  The  FEL  parameters  used  are  chosen  to  be 
illustrative,  and  do  not  represent  a  specific  FEL  design.  In  particular,  the  number 
of  periods,  N  =  3,  is  chosen  to  limit  the  number  of  "loops"  in  the  picture  so  that  the 
motion  would  be  clear.  The  optical  field,  la  I  =  3,  is  below  saturation,  and  large 
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enough  to  cause  noticeable  motion.  The  current,  j  =  1,  was  chosen  to  be  in  the  low 
gain  regime,  so  that  the  model  conditions  would  be  met.  The  motion  demonstrates 
the  same  fast  oscillation  previously  described.  The  point  about  which  the  electron 
oscillates,  however,  evolves  in  a  manner  similar  to  the  the  way  the  electrons  in  the 
helical  imdulator  did  in  the  presence  of  optical  fields.  This  suggests  a  reasonable 
means  for  comparison  of  the  numerical  solution  to  the  approximate  analytical 
result  that  led  to  a  modified  pendulxun  equation  (40).  If  the  center  point  of  the  fast 
oscillations  can  be  followed,  then  the  equation  of  motion  of  the  center  should 
resemble  the  solution  to  the  pendulum  equation.  This  can  provide  a  validation  of 
the  approximate  results,  and  of  the  Bessel  function  coupling  coefficient,  e7o(^>-T/i(^). 


Figure  7.  Exact  numerical  solution  to  motion  with  light  present. 

The  center  of  the  fast  oscillations  can  be  determined  by  the  time-averaged 
phase-space  coordinates  over  each  period  of  the  fast  motion.  The  fast  motion  is 
caused  by  imdulator  deflection,  so  that  the  fast  oscillation  has  a  period  of  exactly 
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half  the  undulator  period.  This  period  can  be  found  exactly,  becaiise  the  average 
longitudinal  velocity  vdthout  light  is  known  exactly,  either  &om  (26)  or  determined 
numerically.  Figure  8  shows  the  results  of  half-period  averaging  of  the  exact 
phase-space  motion  with  lif^t  present.  As  can  be  seen  easily,  the  averaged  motion 
follows  the  exact  evolution.  The  fast  motion  doesn't  need  to  be  followed  within 
each  period,  because  the  optical  fields  evolve  slowly  over  several  periods.  The 
phase-space  trcgectories  do  not  depend  on  the  fast  motion,  but  the  averaged, 
slowly-evolving  motion.  Therefore  the  fast  motion  is  less-interesting,  but  has  the 
undesired  property  that  it  covers  a  much  greater  scale  than  the  slow  motion. 
There  is  no  loss  of  generality  in  following,  but  omitting  to  display,  the  fast  motion. 
The  optical  fields  depend  only  on  the  slow  motion  which  can  be  put  into  the  wave 
equation  (9).  Figure  9  shows  the  same  situation  as  Figiure  8,  but  omits  the  fast 
display. 


Figure  8.  Exact  fast  motion  and  slow  time-averaged  motion. 
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Figure  9.  Slow  time*averaged  motion  with  fast  oscillation  omitted. 


2,  Comparison  with  the  Modified  Pendulum  Equation 

It  is  now  possible  to  use  this  method  to  solve  for  the  evolution  of  a 
uniform  distribution  sampling  of  electrons  in  phase,  and  use  the  results  to  follow 
their  evolution  in  phase-space.  Figure  10  shows  the  case  with  twenty  sample 
electrons,  which  are  equally  spaced  in  time  as  they  enter  the  imdulator,  and  are  at 
the  exact  resonance  energy.  As  it  appears,  they  follow  trajectories  that  could  be 
described  by  solving  the  penduliun  equation. 

This  can  be  compared  directly  to  the  numerical  solution  of  the  pendiilum 
equation  in  Figmre  11  directly  below  it.  The  motion  has  the  same  general  features, 
but  does  not  have  the  same  amplitude.  The  Bessel  fxmction  coupling  coefficient  can 
account  for  the  difference.  For  this  case,  K  =  1,  the  argument  of  the  Bessel 
functions,  ^  *  iiL*/2(l+/ir*)  =  1/4.  This  yields  as  factor  i(4)  =  0.86.  This 
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factor  modifies  two  parameters:  the  optical  field,  Og,  and  tot'  current  j.  With  the 
substitutions,  a  -» a[J<,(^Wi(^)]  and  j  -*  the  new  parameters 

become  Og  =  2.6,  and  j  -  0.71.  Figure  12  is  the  numerical  solution  of  the 
pendulum  equation  using  the  Bessel  function  coupling  coefficient  modifications. 
The  motion  is  now  virtually  indistinguishable  firom  the  exact  solution  in  Figure  10. 
For  these  parameters,  the  Bessel  fimction  coupling  coefficient  modification  to  the 
pendulum  equation  is  an  acciuate  modol  to  the  exact  time-averaged  phase-space 
motion  in  the  linearly-polarized  undulator. 
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Figure  10.  Exact  time-average  numerical  solution  for  linear  undiilator 
with  electrons  exactly  on  resonance,  and  uniform  longitudinal  distribu¬ 
tion. 
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Figure  11.  Numerical  solution  to  the  pendulum  equation  with  identical 

FEL  parameters. 

3.  Range  of  Validity  for  Jo(.%hJii%) 

Although  the  success  of  modeling  the  time-averaged  exact  motion  with 
the  Bessel  function  coupling  coefficient  modification  to  the  pendulum  equation  is 
encouraging,  it  has  only  been  demonstrated  in  a  regime  where  Kly  is  small.  This, 
however,  was  precisely  where  the  approximate  analytical  solution  could  be  found. 
Now  the  question  arises,  what  is  the  range  of  validity  for  the  approximate  result? 
Previously,  there  have  been  few  means  available  to  answer  this  question.  The 
numerical  solution  of  (33)  with  time-averaging  provides  a  visual  means  for 
comparison.  Just  like  what  was  done  previously,  different  values  of  K/y  can  be 
iised,  and  compared  to  the  numerical  solution  of  the  Bessel  function  coupling 
coefficient  modified  pendiiliun  equation.  If  the  motion  is  identical,  then  the  results 
are  compatible. 
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Figure  12.  Numerical  solution  to  the  pendulum  equation  with  Bessel 
function  coupling  coefficient  modifications  to  and  j. 

However,  K/y  cannot  be  chosen  to  be  arbitrarily  large.  The  physical 
significance  of  Kfy  can  be  seen  from  the  previous  results  with  no  li^t  present.  The 
larger  the  magnetic  field,  the  larger  is  the  transverse  deflection  of  the  beam  in  the 
undulator.  This  in  turn  reduces  the  longitudinal  component  of  velocity.  At  the 
most  extreme,  the  deflection  can  completely  turn  the  beam  around,  and  no 
electrons  will  pass  through  the  undulator.  If  y  in  turn  is  made  larger,  the  beam 
will  still  have  enough  energy  to  make  it  past  the  imdulator  deflection.  A  limiting 
case  can  be  found  by  setting  p„  =  l-(l+2K'*yy*  =  0.  For  a  given  value  of  K,  the 
minimum  energy  needed  to  make  it  through  the  undulator  is  Ymi«  =  2^^+l.  For  an 
example  of  IT  «  1,  this  leads  to  and  K/y=  l/Vs  -  0.58. 

At  low  energies  that  are  above  Ymin.  the  resulting  deflection  may  cause 
the  beam  to  transition  from  the  relativistic  regime.  In  phase-space,  this  will  cause 
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some  degree  of  asymmetry.  In  Figure  13,  £*  =  1  and  y  =  2,  and  the  optical  field  is 
Oo  ~  0.  The  resulting  orhits  without  any  li^t  present  are  highly  asymmetric 
because  varies  at  the  top  and  bottom.  This  suggests  that  this  regime  may  be  in 
dear  violation  of  some  of  the  assumptions  used  in  the  standard  derivation  of  (40). 
The  electron’s  velocity  in  the  z  direction  can  not  be  treated  as  a  constant  with 
valuh  ^2  ~  precise  variation  of  oi^t  be  included,  making  it  even  more 

difficult  to  obtain  simple  integral  relations  that  lead  to  the  Bessel  function  coupling 
coeffidents.  An  exact  resiilt  is  impossible,  not  only  in  practice,  but  for  reasons 
discussed  earlier,  the  method  cannot  be  considered  strictly  rigorous. 


Figure  13.  Phase-space  orbits  without  li^t  in  large  if /y  regime. 


Given  that  the  analytical  difficulties  cannot  be  resolved,  the  range  of 
validity  for  the  Bessel  function  coupling  coeffident  can  be  evaluated  by  comparison 
with  simulations.  The  region  of  large  Kly  is  the  most  likely  regime  for  discrepancy. 
Figure  14  shows  the  phase-space  evolution  with  liglit  in  the  high  if /y  regime.  Not 
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only  do  the  orbits  take  on  a  different  shape,  but  the  electrons  appear  to  travel 
further  along  their  orbits  than  before.  An  unexpected  benefit  of  such  motion  is 
unproved  gain  compared  to  the  low  A'/y  model  using  the  penduliun  equation.  Since 
the  pendulum  equation  model  has  no  dependence  on  if/y,  it  is  immediately 
apparent  that  the  two  models  disagree.  The  disagreement  cannot  be  resolved  by 
replacing  the  Bessel  function  coupling  coefficient  with  some  more  accurate 
functional  form,  which  was  the  original  aim  of  this  research.  A  more  fundamental 
problem  lies  in  the  pendviliun  equation  description. 


Figure  14.  Exact  averaged  phase-space  trajectories  for  the  linearly- 
polarized  imdulator  in  the  large  Kfy  regime. 

In  order  to  fully  understand  what  happens,  it  is  necessary  to  return  to 
the  helical  undulator.  Theoretically,  this  model  has  none  of  the  complications  that 
plague  the  linear  case.  The  helical  case  can  be  solved  exactly  with  only  the 
assumption  of  low  gain,  just  as  (33)  was  derived.  The  result  is 


(41) 


Pz  _  AK  [(Dp  -HoC  1-p^  )]cos[(^  -^-kp  )z  -cat  +([>] 

1-pz*  ”  UK^+A^-2AK8m[(k+ko)z-<i)t-¥^] 

This  result  can  be  solved  numerically,  transformed  directly  into  phase-space 
variables  and  displayed  as  in  the  linear  case.  The  complication  of  averaging  is 
unnecessary.  Figure  15  shows  a  large  K/y  case  in  a  helical  undulator.  This,  too, 
does  not  agree  with  the  pendiilum  equation,  shown  in  Figure  12,  which  is 
independent  of  y.  The  assumptions,  K/y«l  and  y=  constant,  used  to  derive  the 
pendulum  equation  do  not  apply  in  this  regime,  and  the  pendulum  equation  is  not 
the  appropriate  description. 


-ii/2  ^  3n/2  0  T  1 


Figure  15.  Phase-space  trsgectories  in  a  helical  undulator  in  the  large 
^/y  regime. 

The  result  that  the  pendulum  equation  itself  doesn’t  apply  in  the  large 
K/y  regime  supports  the  contention  that  an  analytical  result  for  all  regimes  is  an 
impossible  goal.  The  validity  of  the  Bessel  function  coupling  coefficient  doesn’t 


37 


extend  into  this  regime,  not  because  it  is  an  inaccurate  functional  form,  but 
because  there  is  no  appropriate  coupling  coefficient  concept.  That  would  require 
the  same  basic  equation,  but  with  a  different  constant  multipbcative  factor.  The 
basic  equation  itself  is  inadequate. 

If  the  value  of  Kly  is  chosen  such  that  the  pendulum  equation  model  is 
an  accurate  description  for  the  helical  case,  the  accuracy  of  the  linear  model  with 
the  Bessel  function  coupling  coefficient  can  be  evaluated.  By  exploration  with 
different  parameters,  Kly-VA  is  the  limit  of  validity  for  the  pendulum  equation  in 
the  helical  case.  With  identical  parameters,  the  Unear  model  is  compared  to  the 
pendulum  equation  with  the  Bessel  function  coupling  coefficient.  In  Figure  16,  the 
result  for  the  limiting  case  can  be  seen.  The  Bessel  function  coupling  coefficient 
accurately  compares  with  the  exact  solution.  This  leads  to  a  surprising  conclusion: 
The  Bessel  function  coupUng  coefficient  is  accurate  in  all  regimes  of  if/y  where  the 
pendulum  equation  itself  can  be  said  to  apply.  The  conditions  required  to  make 
the  assumptions  invaUd,  also  change  the  basic  physics  involved.  If  a  restriction  is 
made  to  regimes  where  the  pendulum  equation  can  be  used,  then  the  Bessel 
function  coupling  coefficient  compares  extremely  well  with  the  exact  solution. 
Outside  this  regime,  then  a  more  fundamental  equation  like  (33)  or  (41)  must  be 
used  instead  of  the  pendulum  equation. 
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Figure  16.  Exact  averaged  phase-space  trfgectoiies  in  the  linearly- 
polarized  undulator  at  the  limit  of  validity  for  the  pendiilum  equation. 
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IV.  EVALUATING  FEL  PERFORMANCE  FROM 
ACCELERATOR  SIMULATIONS 

A.  ACCELERATOR  SIMULATIONS 

While  accelerator  technology  and  simulation  is  highly  advanced,  integration  of 
thia  technology  with  FEL  technology  has  often  lagged  behind.  This  is  partly 
because  there  is  no  concise  means  of  describing  what  the  effect  of  beam  design  will 
be  on  a  specific  FEL  design.  A  recent  trend  has  been  to  develop  integrated 
computer  models  that  incorporate  both  accelerator  simulation  and  FEL 
simulations.  The  INEX  code  at  Los  Alamos  National  Lab  (LANL)  was  designed  to 
integrate  PARMELA  accelerator  and  FELEX  free  electron  laser  codes.  This 
accomplishment  can  accurately  model  the  entire  FEL  system  from  beam  creation  at 
the  cathode  to  the  laser  output.  The  price  of  this  achievement  has  been 
complexity.  Both  PARMELA  and  FELEX  require  a  vast  amount  of  input  to  specify 
all  the  components.  The  INEX  code  is  difficult  to  run,  takes  a  long  time  on  a 
CRAY  computer,  and  produces  voliuninous  output  This  is  painstaking  and  costly 
which  makes  it  of  value  only  in  the  characterization  of  mature  designs.  In  practice, 
accelerator  and  FEL  designs  are  pursued  independently.  Beam  quality 
requirements  are  usually  specified  only  vaguely,  and  are  often  at  the  limit  of 
current  technology. 

A  simple  code  is  developed  for  use  with  the  accelerator  simiilation.  It  allows 
quick,  and  easy  interpretation  of  the  beam’s  performance  as  an  FEL  There  is  a 
sacrifice  of  accuracy,  but  it  is  of  improved  use  as  a  tool  for  the  designer,  who  can 
test  an  accelerator  design  on  specific  FEL  designs.  This  is  an  improvement  over 
vague  beam  quality  specifications  which  are  independent  of  specific  FEL 
parameters. 
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B.  BEAM  QUAUTY  EFFECTS  ON  UNDULATOR  TRAJECTORIES 

As  was  seen  in  the  preceding  chapter,  the  motion  of  individual  electrons 
inside  the  iindulator  can  be  solved  by  approximate  analytical  and  exact  numerical 
means.  The  exact  tngectory  of  an  individual  electron  can  be  determined  by  its 
initial  conditions.  This  leads  to  the  important  premise  that  beam  quality  manifests 
itself  as  a  spread  in  initial  conditions.  The  final  exit  conditions  from  the 
accelerator  simulation  are  the  initial  conditions  for  the  FEL  undulator.  Therefore, 
the  six  phase-space  coordinates  of  each  of  the  sample  electrons  provides  the  data 
base  for  beam  quality  predictions.  For  Ng  electrons,  there  will  be  6Ng  coordinates. 
Typical  simulations  use  up  to  20,000  sample  electrons,  and  deal  with  120,000 
coordinates. 

An  accepted  measure  of  beam  quality  from  the  accelerator  community  is 
emittance.  It  has  several  definitions  and  has  been  the  source  of  confusion  in  many 
instances.  Emittance,  as  it  will  be  used  throughout  this  chapter,  is  defined  as  the 
transverse  rms  emittance  [Ref.  24] 

Erms  =  4[^<X^XX'^  -  <XX'>^J  (42) 

where  x  is  a  coordinate,  and  x'  =  dxidz.  The  <  >  denotes  an  average  value  over 
the  sample  electrons.  For  extremely  relativistic  beams,  x'  corresponds  to  an  angle 
with  the  longitudinal  axis,  transverse  emittance  has  two 

components,  and  Sy,  each  defined  by  (42).  The  combined  transverse  emittance 
can  be  foimd  as  tmtu  =  High  energies  make  this  number  lower  due  to 

relativistic  effects,  so  a  normalized  emittance  can  be  defined  by  Snorm  =  This 
will  be  explicitly  stated  when  used.  Another  characterization  of  emittance  is  the 
area  enclosed  by  an  ellipse  in  phase-space  that  contains  a  certain  percentage  of  the 
sample  electrons.  An  example  is  the  "90%"  emittance,  ego%  [Ref.  25],  which  is  the 
area  of  the  ellipse  which  encloses  90%  of  the  sample  electrons.  In  practice,  the  90% 
emittance  is  easier  to  measure  than  erm»-  They  are  related  but  depend  on  the 
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specifics  of  the  distributions  to  be  quantitatively  related  [Ref.  26]. 

The  effect  of  initial  conditions  on  trajectories  in  the  undulator  can  be  found  by 
approximate  analytical  results.  First,  introduce  the  additional  complication  of 
non-ideal  fields.  This  is  included  because  it  is  a  common  feature,  especially  in 
linearly-polarized  imdulators.  The  magnetic  field  of  the  undulator  can  be 
represented  by 

IT  =  R  [  0,  sinikoZ  )cosh(Aoy ),  cosfA^oZ  )sinh(A!oy )]  ,  (43) 

where  =  2it/Xo  is  the  \andulator  wave  number,  and  is  the  undulator  period 
[13].  This  assumption  correctly  accoimts  for  off-axis  fields,  with  their  increasing 
magnitude  towards  the  magnetic  pole  pieces.  This  is  also  more  realistic  than  (17), 
because  it  is  a  solution  to  the  Maxwell  equations  for  a  static  field  whereas  (17)  was 
not.  Sometimes  the  pole  pieces  are  shaped  to  provide  a  similar  field  structure  in 
the  other  transverse  direction  because  of  the  desirable  focusing  properties  [Ref  14]. 
This  is  called  two-plane  focusing,  but  will  be  ignored  here.  It  is  a  minor  matter, 
however,  to  add  it  to  the  resulting  code. 

The  motion  of  the  electrons  can  be  found  in  a  similar  manner  as  the  previous 
chapter.  The  effect  of  beam  quality  can  be  found  by  examining  the  electron 
trajectories  without  significant  effect  from  light.  The  wavelength  of  light  will  be 
assumed  to  be  defined  by  the  resonance  condition,  so  that  the  motion  can  be 
discussed  with  respect  to  the  meaningful  phase-space  coordinates  introduced 
earlier,  (C>v).  The  Lorentz  force  equation  gives  the  following  component  equations 

-eB 

p*  = - 5-[  cosikoZ  )y  sinhikoy )  -  z  siaiko  IcoshfA^y )]  ,  (44a) 

ymc^ 

pR 

Py  = - 5-xcos(AoZ)sinh(*oy )  ,  and  (44b) 

ymc^ 

= -^^isin(AoZ)cosh(^oy)  .  (44c) 

ymc^ 

The  X  component  is  a  perfect  derivative  and  can  be  integrated  immediately 
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fix  (y  »2  )  =  Px  (0)  +  V2— cosh(*oy  )[1-cos(A5oZ  )]  , 

7 

where  K  =  eE'Kgl2Tanc^,  and  Px(0)  is  the  fi^  at  z  =  0.  For  nearly  perfect  injection, 
there  is  little  correction  firom  treating  Px(0)  =  Therefore,  the  x  component 

can  be  written  as 

fix  =-'^-yco8h(Aoy)co8(^o2)  •  (45) 

Substituting  (45)  into  (44b)  gives 

P„  =  -2cko  ^cosHko  z  )8inh(2Aoy ) 

r 

The  fast  oscillating  motion  in  the  oosHJkoZ )  term  is  uninteresting  so  it  is  more  useful 
to  take  the  average  motion  over  several  periods  of  the  undulator.  The  averaged 
equation  is 

K  =  -^^^sinhi2koy)  . 

For  small  deflections  firom  the  z-axis,  use  the  small  angle  approximation 
sinh(2Aoy )  ~  2koy ,  so  that 


c^kX 


This  describes  simple  harmonic  motion  in  the  y-z  plane  which  is  called  betatron 
motion.  Using  the  variable  x  =  cf /L ,  this  can  be  rewritten  in  dimensionless  form 

OO 

y  (x)  =  -o)|y(x)  ,  (46) 

where  the  open  circles  denote  derivative  with  respect  to  x,  and  Op  =  2nNKly  is  the 
dimensionless  betatron  frequency. 

The  exact  betatron  motion  is  determined  by  the  initial  conditions.  The 
solution  for  the  y-component  is 
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(47) 


>'(t)  =  >oCos(a)px)  + 


LOy 

COp 


sin((OpT) 


where  is  the  initial  y  position,  L  is  the  undulator  length,  and  0y  is  the  initial 
angle  with  respect  to  the  z  axis  [Ref.  17].  The  betatron  motion  causes  a  periodic 
deflection  in  the  undulator  similar  to  the  periodic  deflection  in  the  x  direction  but 
with  much  longer  period.  Because  the  motion  is  coupled  by  conservation  of  energy 
through  Y”^  =  1  “  Px  ~  Py  “  Pz^j  the  deflection  in  the  y  direction  causes  changes  in 
Pa .  Changes  in  the  x  position  have  no  effect  without  two-plane  focusing.  A  change 
in  the  angle  6^  adds  a  constant  value  to  P^.  The  approximate  average  z  motion  is 

f.  =  1  -  -  fe"  .  MS) 


where  6^  =  .  The  effect  on  v  can  be  computed  from  the  definition 

V  =  L[{.k+ko)^~k].  If  this  definition  is  taken  as  a  transformation  then  for  each 
sample  electron  there  will  be  an  associated  phase  velocity  given  by 


Vj  =L(.k+ko) 


2f 


9/ 

2“P  ,2  “  2 


{k+ko) 


(49) 


where  j  denotes  the  coordinates  of  the  electron.  The  real  resonance  condition 
is  defined  by 


1 

=  if2:v,=o 

0 


(50) 


This  can  be  satisfied  by  a  unique  value  of  k  which  corresponds  to  the  resonant 
optical  wavelength  for  this  electron  beam. 

Now,  a  complete,  dimensionless  description  can  be  made  which  expresses  the 
effect  of  initial  conditions  on  the  phase  velocity  using  (49).  A  sampling  of  electrons 
can  be  used  to  make  a  distribution  plot  of  phase  velocity,  which  in  turn  can  be  used 
as  a  diagnostic  for  PEL  performance. 
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C.  PHASE  VELOCITY  DISTRIBUTIONS  FROM  INITIAL  CONDITIONS 


Using  equation  (49)  as  a  formula  for  computing  the  effect  of  initial  conditions 
on  the  phase  velocity,  each  electron  can  be  assigned  a  value,  Vj ,  as  a  function  of  its 
six  phase  space  coordinates.  This  can  be  made  into  a  distribution  for  each  slippage 
distance,  the  distance  that  light  moves  ahead  of  the  electrons,  since  that  is  the 
relevant  interaction  length  in  the  micropulse.  However,  there  may  be  thousands  of 
slippage  lengths  in  a  long  micropulse  with  a  short  optical  wavelength,  thereby 
making  the  number  of  plots  prohibitive.  A  short  characterization  of  each  slippage 
length  may  be  given  by  the  centroid  and  rms  width  of  each  distribution.  If  the 
basic  shape  of  the  individual  distributions  are  roughly  the  same,  then  this  will  be  a 
useful  simplification. 

Figure  17  shows  a  t3rpical  slippage-length  phase-velocity  distribution.  The 
resonance  parameter  refers  to  v.  Since  exact  resonance  is  defined  by  v  =  0,  this  is 
effectively  a  plot  of  Av.  The  sample  electrons  came  from  the  PARMELA  simulation 
of  an  accelerator  design  for  a  compact  FEL  imder  development  at  LANL  Beam 
parameters  are:  peak  current,  1  »  400A;  beam  energy,  =  ISMeV;  and  pulse 
length,  =  16ps.  The  FEL  parameters  are:  undulator  period,  Xo  =  0.9cm;  number 
of  periods,  N  =  15;  and  imdxilator  parameter,  K  =  0.7.  The  resonance  condition, 
(50),  gives  an  optical  wavelength  of  0.7pm.  Beam  quality  was  on  the  average 
e  s  lOx  mm-mrad,  and  energy  spread  0.07%. 

The  shape  of  Figure  17  is  highly  asymmetric.  This  can  be  understood  in 
terms  of  dividing  the  effects  into  two  broad  categories:  transverse  position  and 
angle,  and  real  energy  spread.  Any  combiimtion  of  initial  transverse  phase  space 
coordinates  will  translate  into  a  reduction  in  the  phase  velocity.  Specifically,  the 
resulting  betatron  oscillation  constitutes  a  reduction  in  the  parallel  component  of 
velocity,  and  therefore  results  in  a  reduction  of  the  phase  velocity.  Real  energy 
spread  may,  or  may  not,  have  a  well  defined  shape.  Most  typically,  whatever 
effects  that  result  from  unequal  acceleration  tend  to  be  symmetric.  Therefore,  a 
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Figure  17.  Phase  velocity  distribution  in  one  slippage  distance. 

beam  dominated  by  real  energy  spread  would  have  a  more,  or  less  symmetric 
distribution  while  a  beam  dominated  by  angular  divergence  would  show 
pronoimced  inhomogeneous  broadening.  An  approximation  to  the  distribution  dut 
solely  to  angular  divergence  is  a  straight  exponential  The  beam  in  Figure  17 
is  an  example  of  a  case  dominated  by  angular  divergence.  For  the  beam  used,  the 
real  energy  spread  is  AyYsO.07%,  and  the  transverse  rms  emittance  is 
es  lOx  mm-mrad.  Without  specific  FEL  parameters,  these  quantities  alone  cannot 
be  interpreted  usefully. 

Figure  18  shows  the  rms  width  of  the  distribution  in  each  slippage  distance 
as  a  function  of  the  position  along  the  micropulse.  In  the  center  of  the  micropxilse, 
the  value  remains  essentially  steady  at  Av  s  2.0.  This  value  must  be  compared  to 
the  relevant  small-signal  gain  bandwidth,  Av^  x,  as  shown  in  Figure  5.  This 
distribution  would  fall  entirely  within  the  gain  bandwidth,  and  therefore,  one 
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would  not  expect  any  serioxis  gain  degradation.  This  is  not  the  only  comparison 
that  can  be  made.  It  is  also  possible  to  take  the  resulting  distributions  and  use 
them  with  the  theories  of  beam  quality  effects.  Namely,  it  is  a  simple  extension  to 
take  the  res\ilts  and  apply  them  to  the  theory  of  Blau  and  Colson  and  find  a 
characteristic  function  that  can  be  used  in  the  "integral  equation"  [Ref  27].  Since 
this  area  has  already  been  well  explored,  it  is  useful  to  take  advantage  of  the 
previous  work  as  much  as  possible.  The  phase  velocity  distributions  provide  the 
bridge  to  this  useful  area. 


Figure  18.  Phase  velocity  distribution  rms  variation  as  a  function  of  position 
along  the  micropulse. 
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The  shape  of  figiire  18  also  shows  a  dramatic  drop  in  beam  quality  at  the  end 
of  the  micropulse.  However,  trends  in  this  data  must  be  compared  to  the  particle 
density.  Large  variations  at  the  ends  are  mitigated  by  the  relatively  low  currents 
there.  Figure  19  shows  how  particle  density  varies  along  the  micropulse.  The  ends 
have  relatively  low  current,  and  therefore,  temper  the  effects  of  beam  quality  there. 


The  way  phase-space  coordinates  are  manipulated  by  accelerator  beam-line 
components  and  focusing  elements  can  dramatically  affect  PEL  performance  even 
though  the  inherent  beam  quality  remains  the  same.  Liouville’s  theorem 
guarantees  that  the  six-dimensional  phase-space  hypervoliune  will  remain  constant 
for  conservative  forces  [Refs.  21,25].  Because  the  beam  is  highly  relativistic,  only 
the  transverse  phase-space  distribution  can  change.  Liouville’s  theorem  then 
predicts  that  the  transverse  emittance  is  a  constant  of  motion.  Although  the 
emittance  is  constant,  different  combinations  of  the  phase-space  coordinates  do  not 
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guarantee  the  same  effect  in  the  FEL.  It  has  been  shown  that  transverse  phase- 
space  coordinates  must  be  "matched"  to  achieve  the  best  performance  for  a  given 
inherent  beam  quality.  The  matching  condition  [Ref.  26]  is 


<y^  = 


eL 

27C0)p 


(51) 


While  this  can  be  done  for  any  small  fraction  of  the  micropulse,  it  cannot  be 
achieved  for  the  entire  beam  simultaneously.  Some  middle  ground  must  be  chosen 
to  hopefully  achieve  optimvim  FEL  performance. 

This  also  points  out  the  inacciiracy  of  using  a  single  transverse  emittance 
value  for  the  entire  pulse.  If  we  look  at  the  transverse  emittance  at  any  one 
slippage  length,  it  remains  at  a  relatively  steady  value,  consistent  with  liouville’s 
theorem.  However,  the  actual  phase  space  plots  may  be  distributed  in  radically 
different  ways.  Although  the  plots  take  up  the  same  area,  an  ellipse  enclosing  a 
msqority  of  the  points  can  be  highly  eccentric.  Since  emittance  measiires  only  the 
area,  a  matched  ellipse  by  (51)  has  the  same  emittance  as  a  highly  eccentric 
ellipse,  although  the  later  doesn’t  satisfy  (51).  The  unmatched  section  may 
oscillate  on  the  whole,  whereas  the  matched  section  is  balanced  and  remains  in  the 
same  distribution  throughout.  Transverse  emittance  as  a  function  of  longitudinal 
position  in  the  micropiilse  is  shown  in  Figure  20,  and  provides  no  information 
about  how  the  individual  sections  are  matched.  If  all  the  slippage  lengths  are 
projected  onto  the  same  plot,  the  resulting  emittance  may  be  much  greater  than 
any  individual  slice.  This  is  because  there  may  be  large  areas  of  non-overlapping 
coverage  that  together  enclose  a  larger  area  than  any  one  slice.  A  single  number 
that  measures  the  transverse  emittance  of  the  beam  cannot  account  for  any  of 
these  effects,  nor  does  the  plot  of  emittance  as  a  fimction  of  position  in  the 
micropulse.  The  only  way  to  incorporate  all  of  these  considerations  is  to  examine 
the  phase  velocity  distribution  functions  at  each  slice,  or  more  simply  a  plot  of  the 
rms  width  of  these  distributions  as  a  function  of  position. 
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Figure  20.  Transverse  emittanoe  as  a  function  of  position  in  the  micropulse. 

In  addition  to  the  nns  width  of  the  distributions,  the  maximiun  or  centroid  of 
the  distribution  may  vary  from  slippage  distance  to  slippage  distance.  Figure  21 
illustrates  the  result  obtained  finom  the  accelerator  design  under  review.  Again, 
there  is  a  strong  variation  at  the  ends.  The  same  argument  presented  previously 
applies  to  the  extreme  variation  at  the  ends.  However,  there  may  be  a  systematic 
variation  along  the  entire  length  of  the  micropulse.  In  systems  where  the  slippage 
length  is  short  compared  to  the  micropulse,  this  may  manifest  itself  as  a  temporal 
firequency  variation  in  the  optical  pulse.  However,  this  effect  is  not  limited  to  a 
frequency  "chirp,"  but  the  details  of  how  this  works  in  short  pulse,  long  slippage 
distance  systems  which  account  for  collective  effects  has  not  yet  been  addressed. 
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Figure  21.  Variation  of  phase  velocity  along  micropulse. 


In  summary,  the  six  phase-space  coordinates  of  all  the  sample  electrons  can 
be  put  into  a  noxninal  nxunber  of  plots  that  describe  much  of  the  essential  FEL 
physics.  The  particle  density  show  what  weight  to  give  the  various  effects.  The 
average  phase  velocity  describes  any  systematic  variation  in  the  frequency  in  the 
optical  pulse,  but  slippage  length  communication  must  be  considered.  The  width  of 
the  phase  velocity  distribution  along  the  micropulse  most  clearly  illustrates  the 
beam  qualify  effects,  and  can  be  used  with  any  number  of  well  established  theories. 
Various  standard  quantities  can  be  computed  for  comparison  to  the  more 
frequently  available  measurements,  for  example,  those  obtained  by  experiment. 
And  finally,  this  program  can  easily  recommend  values  for  beam  matching 
parameters. 
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